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This paper shows that the backscattering of electromag- 
netic and gravitational waves can be dominant when the radi- 
ation is produced very close to a spherical black hole. Numeri- 
cal investigation shows that almost 50 percent of the outgoing 
quadrupole gravitational wave is backscattered, for a class of 
initial data. A similar analysis reveals at least 20 percent ef- 
fect for a dipole electromagnetic radiation. Numerical results 
confirm theoretical predictions that the backscatter of short 
wavelength radiation is negligible. In the long-radiation band 
is observed a rather weak dependence on the wavelength. Our 
studies base on the linear approximation. They can be of rel- 
evance for the determination of the total energy of backscat- 
tering tails and quasinormal modes. 
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I. INTRODUCTION 

This paper is dedicated mainly to a numerical inves- 
tigation of the backscattering ( [Q - (TIJ) of waves in a 
Schwarzschild spacetime. It supplements the earlier an- 
alytic studies of scalar Eia, electromagnetic ( |lif| , 14 1) 
and gravitational fields j[16[). We study the propagation 
of the electromagnetic and ( even-parity) gravitational 
waves in a background Schwarzschild spacetime. It is as- 
sumed that initial data describe an isolated pulse (burst) 
of a (gravitational or electromagnetic) wave. As in the 
earlier studies ( [fL2| , |l3} | and 16 1) the strength of the 
backscattering is assessed by finding the fraction of the 
initial burst energy that will not reach a distant observer 
in the main pulse. The numerical investigation, that is 
reported below, gives a quantitative evaluation of the ef- 
fect. 

Sec. II brings an analytic estimate of the backscatter- 
ing of the electromagnetic fields. The obtained result is 
a substantial improvement of a former bound p4] |. Sec. 
Ill is dedicated to the presentation of numerical results 
on the backscattering of electromagnetic fields. The ef- 
fect depends both on the relative width of the initial data 
and on the distance. Sec. IV describes the propagation of 
quadrupole gravitational waves in a Schwarzschild geom- 
etry. Relevant energy formulae are defined and analytic 
estimates are reminded. Section V goes on with numeri- 
cal studies of the backscatter of gravitational fields. Sim- 
ilarly as before, the backscattering happens to depend on 
the relative width of the initial data and on the distance. 



In Sec. VI we briefly review the obtained results. 

The space-time geometry is defined by 
Schwarzschildean line element, 
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where t is a time coordinate, R is a radial coordinate that 
coincides with the areal radius and d£l 2 = d9 2 + sin 2 0d<p 2 
is the line element on the unit sphere, < <p < 2ix and 
< 8 < ir. Throughout this paper G, the Newtonian 
gravitational constant, and c, the velocity of light are 
put equal to 1. We define the Regge- Wheeler coordinate 
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i? + 2mln( 
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and r\ r = 1 — 2m/r. 



II. BACKSCATTER OF ELECTROMAGNETIC 
WAVES: NEW ANALYTIC ESTIMATES 

We seek, following (Of, a solution <3>(r* ,t) in the form 



* = * + S, 
where S satisfies the dipole equation 



(-d 2 + d 2 ,)s = VR 



R 2 



Here 



*(r*,i)=fi t /(r*-i) + 
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R 



(3) 



(4) 



(5) 



and / is an arbitrary function with support in (a, oo). / 
can be uniquely determined from initial data correspond- 
ing to an initially outgoing radiation. * solves Maxwell 
equations in Minkowski spacetime and it corresponds to 
the dipole radiation. Initially 6 = OqS = 0. 

The energy En(t) of the electromagnetic field ^ con- 
tained in the exterior of a sphere of the radius R reads 



E R (t) 
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E a = E a (0) is the energy of the initial pulse. Let an 
outgoing null geodesic T r originate from a point (r, 0) of 
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the initial hypersurface. In the Minkowski spacetime the 
outgoing radiation contained outside T a does not leak 
inward and its energy remains constant. In a curved 
spacetime some energy will be lost from the main stream 
due to the diffusion of the radiation through T a . 

The energy loss is equal to a line integral along r a 
(where / = # = 0), 



The integrand of the external integral with re- 
spect "t" is equal to the directional derivative 



1/2 



and the last line of the pre- 



ceding equation is equal to 
d 
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An energy H (R, t) of field 6, 

H(R,t)= (°° dr(^^ + Vr (d r Sr + (S) 2 ^), (8) 
Jr v Vr r z J 

satisfies the "conservation" equation ju| 
(d t + d* r )H(R,t) = 
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The integration of (|^) along f a yields 

#(oo,oo) - H(a,0) = - 12m / dt / drd a 5^-; 

27r Jo J a t r 

(10) 

here (a t ,t) S T a and _ff(oo,oo) is the asymptotic energy 
of the field 6. Eq. (|l0|) implies (taking into account that 
H(a, 0) = 0, due to the vanishing of initial data of the 
field 6) that 
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This integral is bounded from above by 



dty/H(t) / dr 
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On the other hand, Eq. (m implies 
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Therefore, (0|15j) imply 



SE a < 36m 2 



dt 



dr 
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One can show, using reasoning as in jlj], that on the 
initial hypersurface t = 



J(R,t = Q) 
1 i?-2m 1 



r — 2m 



dr 



(r — 2m) 2 
(17) 



this can be written, using the Schwarz inequality, as 
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Notice that J,f dr(f ) 2 < S a /(4?r). From (jl|) follows 



< E a amF(fh,y), 



(19) 



where 

F(m, y) = y — 1 + 



16m 



16m 
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Here m = m/a and ?/ = R/a. Inequality ( ppf ) constitutes 
a refined version of Eq. (13) in flTij . 

Let (i?,i) S f r . The right hand side of (§(]) is an 
increasing function of r; since R > r and / is constant 
along r r , one obtains 



4^/ 2 (i?,i) 



R2 ' - E a ai] R F(m,R/a), 



(21) 
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Define x = a t /a. The integrand of the right hand side of 
( |l6| ) is bounded, taking into account (pl|): 



2tt 



E a a 
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dr = ^±G(m,x). 
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Here G(m,x) — J°° dx(x — 2m) 3 F(m, x)/x 9 ; the rather 
lengthy integration can be done explicitly with the use 
of Mapple. The result is too long to be presented here.. 
The insertion of (^|]) into ( |l6| ) leaves us with 



SE n < 
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dx 
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Here we replaced the time variable t by the radial variable 
r; notice that along T a one has dt — drjr\ r . The use of 
the Schwarz inequality allows one to bound J23|) from 
above by 



8E a < 



l%m 2 E a 



dxG(rh, x)x 2 



dx 
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That implies the following estimate for the ratio of the 
diffused energy by the initial energy of the pulse of radi- 
ation 



SE a 
E a 



where 
C(m) 
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C(m) 
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828m 4 



10080((-l + 2m) 4 m 4 ) 

44m 3 + 2352m 6 + 2016m 4 ln(l - 2m) + 
2688m 6 ln(l - 2m) - 4032m 5 ln(l - 2m) - 
360m 3 ln(l - 2m) + 6m + 36m 2 ln(l - 2m) 
-30m 2 + 31n(l - 2m) - 

18mln(l-2m)^ (26) 

This result can be improved in the case when the initial 
pulse is located outside the sphere a = 10m/3. One finds, 
eliminating the time derivatives of 'J, that the energy can 
be written as 



E a {t)=4irJ dr\^ r (d r ^) 2 + 

29^(0.9-^)1 

r>0 if* I 
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(27) 



Thus J a °° ^ 2 /r 2 dr < £ q /(8tt) if a > 10m/3. As a con- 
sequence the factor in the inequality (^) halves and, fi- 
nally, the formula (|25|) can be replaced by 



5E a /2m\2 , 

-^<2.25(-) Cirn) 



(28) 



A more careful estimate, in which the d r ^-related term 
is taken into account, gives even a stronger result, with 
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(30) 



that significantly improves the former result proven in 
|l4j| . The coefficient C(m) diverges at m = 0.5 but 
it depends rather weakly on m = m/a in the range 
(0,0.25). For instance, for m/a = 0.1 and m/a = 0.25, 
one has 2.25C(m)/(l+0.125?7 a ) w 0.39 and 2.25C(m)/l+ 
0.125?7 a ) w 1.4, respectively. 

It can be of interest to consider the case of initial data 
of compact support (a, 6). The simplest estimate can be 
obtained as follows. It is easy to see that one has 



^f 2 (R,t) 
R 2 



< E a ar]%F(m, b/a), 



(31) 



instead of (^l|). The insertion of (|T|) into Eq. (|T^) yields, 
after a calculation analogous to that performed above, 
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This is valid also inside the photon sphere. In the case 
when (6 — a) /a << ry 3 , one can clearly see that the 
backscatter is negligible. 

A simpler and much stronger estimate can be obtained, 
assuming that 2m << a and the support of initial data 
satisfies the condition b — a << a. In this case one has 
from ((l|) 



(32) 
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Here vE'(a) = and one can show, employing the same 
approach as before, that 

»R P R 

drd r M < 
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Thus 



^f 2 (R) < Ea b2 

R ~ 4 



(i-t) 3 



(34) 



(35) 



The insertion of ( |35| ) into (|T|) yields (notice that the 
spatial integration now extends from a t to bt and that 
(a* - b t )/a ps (6 - a)/a) 

8E a / m \ 2 / 6 — a \ 4 



0£j a /m\' 
— - < 0.51 — ) 

E n . ~ \aJ 



(36) 
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III. NUMERICAL RESULTS: 
ELECTROMAGNETIC FIELDS 

This section reports numerical results on the backscat- 
ter. We begin with demonstrating that the effect can 
be strong and then discuss qualitatively its wavelength 
dependence. 

We choose initial data generated by the function 

^2^ 



dtf 



(R- 
= 



a) 2 - 01 exp( 



-w{R 
2m < R < a. 



R> a: 



(37) 



There are two free parameters, a and w. The exponent 
2.01 guarantees that the energy density (that depends, 
in particular, on d\f) vanishes at a. 

The initial data have a noncompact domain, but they 
are very small outside a compact support. In particu- 
lar, the projected one-dimensional initial energy density 
(the integrand of Eq. (||)) is practically of compact sup- 
port. On the other hand, the function / is obtained from 
Eq. ( |37| ) by integration and it is constant asymptotically. 
Since / enters the evolution equation , it is responsible 
for the backscattering; its asymptotic constancy ensures 
that the effect is relatively strong. In what follows the 
mass m is assumed to be 1. The backscattered energy is 
obtained by a numerical approximation of the line inte- 
gral (§. 

We integrated Eq. (||) for a = 2,4 and 20. Fixing the 
parameter a, calculations were performed for a number 
of w's. The results seem to imply that 8E a /E a has a 
single maximum at some u>*, with SE a /E a being almost 
constant in some vicinity of w» . The limitation inherent 
to the numerics does allow us to find only an approximate 
value of u>* ; that value will be called as a " maximal" point 
and the corresponding value of 6E a /E a will be referred 
to as a " maximum" . That terminology will be kept in 
the rest of Sec. Ill as well as in Sec. V. 



A. The maximal backscattering 

The strongest backscattering effect is found with the 
choice of parameters a = 2 and w = 5 x 10 6 . The 
backscattered energy has been obtained by the integra- 
tion of (0) along r 2 .oooi @- Under these conditions 
the value of the ratio 5E a /E a approximates 20.5%. The 
backscatter is relatively insensitive on the choice of w. 
Fig. la) shows the initial energy density for various val- 
ues of w as well as the strength of the backscatter as 
measured by the ratio 5E a /E a . 



B. Resonant type initial data and relative-width 
dependence 

Figs, lb) and lc) show the behaviour of the initial 
energy density that corresponds to a = 4 and a = 20, 



respectively, and to a number of values of w. A number 
assigned to a curve gives the ratio of the backscattered 
energy. It is obvious that in these numerical examples the 
strength of the backscatter is correlated with the width of 
the support A of initial data. The relative width A/a of 
the energy density is related to its spectral composition; 
that suggests that the effect depends on the wavelength 
of the initial radiation. 

Numerical investigation confirms a theoretically de- 
rived conclusion that in the limit A/a — > there is 
no backscatter ( see also Eq. (|3"rj|)). In our ex- 

amples small values of the ratio A/a correspond to 
w >> 1. We found that the ratio SE a /E a monotoni- 
cally decreases with w, if w is large enough. In order to 
demonstrate how dramatic the changes can be, we com- 
pare data concerning the case a = 4. If w — 10 3 , then 
8E a /E a = 4.5 x 10~ 8 ; that is rougly 10 4 less than for 
w = 1, when SE a /E a = 1.9 x 10~ 4 . 

The values w = 0.1 (in the case of a = 4m) and 
w = 0.001 (in the case of a = 20m) correspond to "res- 
onant" initial data. That "resonant" behaviour is rather 
weak; for instance in the case of a = 4 the change of 
w by a factor of 10 from the maximal case results in a 
relatively small, less than 50%, change of the efficiency 
factor SE a /E a of the backscattering. In general, we can 
say that the "resonant" width of the initial energy density 
(which gives also the order of the " resonant" wavelength 
of the radiation, if a >> 2) is of the order of the distance 
of the radiating source from the black hole. 

Fig. la) clearly shows, in contrast with the preceding 
cases, that the strength of backscatter is rather weakly 
correlated with the width of the support of initial data. 
The reason for that is that (for this set of configurations) , 
the smaller is the width, the smaller is the distance from 
the horizon (and thus the redshift increases). There are 
two competing effects that work each against the other. 

In those cases that have been studied in our paper, the 
backscattering is strongest at a = 2 and w — 5000000, 
when SE 2 . 0001/^2. 0001 = 20.5%. The resonant efficiency 
drops then by a factor of 500 when passing to another res- 
onant case at a = 4 and w = 0.1: SE4/E4 = .045%. This 
dramatic change is expected, since the backscattering is 
strongest within the sphere R = 3m. The largest value 
of 6E20/E20 (and w = 0.001) is obtained for w = 0.001 
and it reads 0.00085%. Thus the increase of a from 4 to 
20 again results in a roughly 50- fold decrease of the reso- 
nant efficiencies, which is somewhat quicker than fall-off 
suggested by the analytic estimate (p5[). 



C. On the comparison of analytic and numerical 
results 



Let us comment on the comparison of the numerical 
data with the analytic estimates of the preceding section. 
The inequality (25, valid generally, is not sharp. In the 
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case of resonant initial data, for instance, it gives a non- 
trivial information only for a = Am and a — 20m (Figs, 
lb and lc); the numerically obtained value of 5E a /E a is 
by three orders smaller than the analytic bounds. The 
formula (|36|), that is valid at large distances and for ini- 
tial data with small relative widths, should be more ef- 
ficient. The use of ([36]) is not allowed, strictly saying, 
in our case since the numerical examples do not satisfy 
the required assumptions concerning the compactness of 
initial data. In the limit of A/a — > 0, however, the initial 
data have " almost" compact support and in this context 
it is interesting that ( |36| ) gives predictions that are com- 
parable with numerical ones, in the case of initial data 
with narrow support. We will study elsewhere the ques- 
tion whether (and under which conditions) the analytic 
criteria are strict. 



IV. PROPAGATION OF GRAVITATIONAL 
WAVES 

We will search the solution of the Zerilli equation fHf ! 
in the form 



* = * 



(38) 



where S is an unknown function satisfying the 
(quadrupole) equation 



R R 2 



2mr)R 
Here 



-3*i 



,£2 
' R 



V(R) = 6 v 2 + m 



63m 2 (l + g ) 

3m \2 



R 2 '"2i? 4 (l + |f) 



(39) 



(40) 



and *i(r* — t), i = 0, 1, 2, are functions that satisfy the 
relations 



ft*i = 3*o 

ft*2 = *i - mftfi. 



(41) 



The combination 



isto(r ._ i)+ M^ + M^i>, (42) 



R 



R 2 



which represents a purely outgoing radiation, solves the 
Zerilli equation in Minkowski space-time (m=0). We 
choose * = *, ft* = ft*, which implies S — ft<5 = at 
t = 0. 

The initial energy density multiplied by R 2 reads 



p = ((ft*) 2 + (ft**) 2 + V^ 2 )/ VR . 



(43) 



The initial data are assumed to be smooth and to be 
nonzero outside a sphere of a radius a > 2m. Thus p is 
smooth and it vanishes on the boundary a. The energy 
content inside a part of a Cauchy hypersurface E t that 
is exterior to a ball of a radius R can be defined as 



E(R,t) 



drp(r, t). 



(44) 



We omit a normalization constant in the definition of the 
energy E(R, i), since we will be interested only in the rel- 
ative efficiency of the backscatter and the normalization 
factor cancels out. The total initial energy corresponding 
to the hitherto defined initial data will be written as E a . 

The energy loss, that is the amount of energy that 
diffused inward T a is equal to a line integral along r o , 



SE a =E(a,0)-E o 



dr 



Vr 



f—(d t + d r ,)S(R,t)) — 



VR 



vs 2 



(45) 



It is necessary to point out that in the case of the ini- 
tial point Rq > a the result would be more complicated; 
the differentiation of the energy along would depend 
also on *o,*i and *2- If, however, the outgoing null 
geodesies is T a , then it starts from a where *o, *i and 
*2 do vanish. Since these functions depend on the dif- 
ference r* — t, their values along outgoing geodesies are 
constant, and that allows one to conclude that they van- 
ish at r a . 

The fraction of the energy that could diffuse through 
the null cone r a satisfies flCf 

Theorem. The efficiency of the backscattering, 
5E a /E a satisfies the inequality 



^<54.5x(^) 2 + 0(m 3 /a 3 ). 
E a V a J 



(46) 



In the case of compact initial pulses one has |l6[ , assum- 
ing (b — a)/a « 1 and m/a « 1, 



E„ 



V a ) 



(47) 



V. NUMERICAL RESULTS: GRAVITATIONAL 
FIELDS 

Below we shall describe numerical results on the 
backscattering of even-parity gravitational waves. As in 
the case of electromagnetic fields, initial data are found 
that give rise to the strongest effect and then the wave- 
length dependence is discussed. Finally, the efficiency co- 
efficient SE a /E a corresponding to resonant initial data is 
shown to decrease with the increase of a distance. 

We choose initial data generated by the function 
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*i = (R - ay- m exp 
= 



-w{R-a) 2 ), R>a; 
2m < R < a. 



(48) 



The exponent is taken to be 2.01, in order to guaran- 
tee that the initial energy density (that involves second 
derivatives of ^i) vanishes at R = a. There are two 
parameters, a and w. The initial data are of noncom- 
pact support but from the numerical point of view they 
are zero outside a compact set. The one-dimensional pro- 
jected initial energy density (^3|) becomes negligibly small 
outside a compact support. The function \&2 is obtained 
by integration (compare ([48]) and (f4l|)) and it is constant 
asymptotically, and \&2 appear in the evolution equa- 
tion (Q). They probably generate the dominant contri- 
bution to the backscattering and the constancy of ^ at 
spatial infinity can ensure the best conditions for having 
the strongest effect. We put m = 1. As in Sec. Ill, Eq. 
(|J) is integrated for a — 2, 4 and 20. The initial energy is 
calculated from the formula (44) and the backscattered 
energy is found numerically from Eq. (45). For each 
fixed parameter a is determined a value w at which the 
factor 5E a /E a is the largest one (but see the explanation 
in Sec. III). 



A. The maximal backscattering 

The highest ratio 8E a /E a is found for parameters 
a = 2 and w = 10 4 . The backscattered energy was cal- 
culated along r2.ooi [0 and then it was noticed that 
SE 2 .ooi /E2. 001 ~ 47%. Let us point out that this exceeds 
by a factor of 10 a prediction made in ]lS| ]. The backscat- 
ter does not depend strongly on w. Fig. 2a) shows the 
initial energy density for various choices of w. A number 
assigned to a particular curve shows the strength of the 
backscatter - the corresponding value of SE2.001/ E 2 .ooi- 



In our numerical examples the ratio A/a becomes smaller 
if w increases. It was found, for sufficiently large values 
of w, that SE a /E a monotonically decreases with w. For 
instance, if a — 4m, then the change of w from 10 _1 to 
4 x 10 3 results in the decrease of 5E a /E a by a factor of 
10 8 : from 3.7 x 10~ 3 to 1.35 x 10 _u . 

The values w = 0.01 (in the case of a = 4m) and 
w = 0.0001 (in the case of a = 20m) correspond to the 
maxima of the factor 5E a /E a . Again, as in sec. Ill, "res- 
onant" peaks are rather mild. In the case of a — 4 the 
change from w — 0.1 to w — 0.001 changes the efficiency 
from 5E 4 /Ei = 0.37% to 5Ei/E 4 = 0.4%, whereas the 
maximum is 0.65% at w = 0.01. Notice, however, that 
the " resonant" width of the initial energy density (which 
gives also the order of the "resonant" length of the ra- 
diation, if a » 2m) is now much bigger (by one order) 
than in the case of the electromagnetic dipole radiation. 

The backscattering is strongest at a = 2 and w = 
10000, when SE 2ml /E 2 . 00 i = 47%. The efficiency de- 
creases quickly, circa 70-fold, when the pulse is moved 
outside the photonic sphere to a = 4. Then the reso- 
nant value 5E4/E4 — .65% corresponds to w = 0.01. At 
a = 20 and w = 0.0001 we found 5E 20 /E 2 o = .01%. 
Thus the increase of a from 4 to 20 results in 65-fold de- 
crease of the efficiency at resonant cases, which is more 
than the fall-off suggested by the analytic bound ([l6|). 

The comparison of the numerical data with the ana- 
lytic estimates ( (46] ) and ( flTf ) yields conclusions similar to 
those made in the electromagnetic case. Again, the width 
independent criterion d46| ) is rather imprecise. In the case 
of a = 20 (Fig. 2c) the numerical value of SE a /E a is circa 
3 orders smaller than that predicted analytically. The 
prediction of the bound (f47j), that is valid for a >> m 
and b — a << a can be more precise. 



VI. CONCLUSIONS 



B. Resonant backscatter, relative width dependence 
and analytic estimates 

Figs. 2b) and 2c) show the initial energy density at 
a = 4 and a = 20, for a selection of values of w, with 
numbers showing values of the fraction of the backscat- 
tered energy that are associated with particular initial 
data. Figs. 2b) and 2c) reveal that the strength of the 
backscatter is correlated with the relative width A/a of 
the support of initial data. In contrast with that, Fig. 
2a) shows a rather weak dependence of SE a /E a with the 
(initial) relative width. The explanation for that anoma- 
lous behaviour is as in the case of Fig. la - that there 
do appear contradictory effects - the diminishing of the 
relative width (with w being increased) goes in pair with 
the decrease of the distance from the horizon. 

The efficiency of the backscatter quickly diminishes 
with the decrease of A/a, as theoretically predicted [M. 



There took place a debate, about decade ago, on the 
quantitative evaluation of the backscattering of gravi- 
tational waves (see Jl8| and references therein). While 
one of the proponents was originally in favour of a very 
strong damping of the radiation, that was implied by the 
backscattering, the final conclusion was that the effect 
is weak and only a small fraction (of the order of a few 
percents) of the long-wave band of the gravitational ra- 
diation can be backscattered Jl8| . In the light of that the 
main results of this paper come as a surprise. It is quite 
likely that the ratio of almost 50% of the backscattered 
gravitational quadrupole radiation, that we find here, can 
be improved for more suitably chosen initial data. The 
same can be said about the backscatter of electromag- 
netic waves, where we established that this effect exceeds 
20 %, for the dipole radiation. 

The numerical examples of Sects. Ill and V confirm 
theoretical predictions that if the relative width of the 
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initial pulse tends to zero then the effect becomes neg- 
ligible. This can be translated (using the so-called sim- 
ilarity theorem of Fourier transforms) into dependence 
on the asymptotic wavelength (or frequency) of the ra- 
diation pj . The simplest argument would invoke to the 
Heisenberg principle, which clearly implies that the com- 
pression of the support of a function leads to the increase 
of the frequency scale in its Fourier transform. The ana- 
lytic proofs can be deduced from formula |6|) and (^) of 
this paper. On the other hand, in the "resonance" regime 
the dependence on the relative width (or, accepting the 
preceding arguments, on the wavelength) is rather weak, 
according to numerical data of Sees. Ill and V. That 
seems to be in a sharp contrast with what is known in 
the case of stationary processes @ . 

Let us point out the difference between the manifes- 
tations of the gravitational redshift and of the backscat- 
tcring. The gravitational redshift is responsible for the 
weakening of the intensity of a radiation, but - barring 
backscatter - all of the initial energy eventually reaches 
an asymptotic observer. In the limit of the geometric 
optics the gravitational redshift is the only phenomenon 
that can be observed. The backscattering in turn is re- 
sponsible for the loss of energy, and it may be important 
in the case of a low-frequency radiation. 

We observed in Sees III and V that the strongest effect 
takes place very close to the horizon. That reflects well 
the known fact that the bulk of the backscatter happens 
inside the marginally stable photon sphere R — 3m. As 
a consequence, the transport of an initial pulse from a 
location close to the horizon to a point behind 3rn must 
be associated with a significant decrease in the efficiency. 
That was in fact observed in our numerical examples, 
both for gravitational and electromagnetic waves. This 
fall-off with a distance can be faster than 1/R 2 , which is 
typical for our analytic estimates. 

The comparison of the analytic and numerical results 
suggests that analytic bounds ((|25|) and (^6|)) are not 
strict; their predictions are bigger from our numerical 
data by two - three orders. That gap can be made much 
smaller by a more suitable choice of initial data, but we 
do not expect that it could be nullified. The estimates 
( |36| ) and ( |47| ) that are specialized to the case of small rel- 
ative widths are expected to be much sharper. The an- 
alytic bounds on the efficiency of the backscatter clearly 
show that the effect is negligible if the initial burst is 
located far from the horizon, irrespective of the detailed 
character of initial data. 

We would like to point out a new important appli- 
cation of the methods developed in this paper. Namely, 
the backscattered (or diffused, in the terminology of |13|) 
energy bounds from above the total energy of tails and 
quasinormal modes of black holes and of a class of quasi- 
normal modes (" w-" modes) of neutron stars ( Ej , ) • 
Therefore the present methods can be used as the first 
step in order to estimate the fraction of the energy of 



gravitational waves that can be carried by quasinormal 
modes in spherically symmetric spacetimes. That would 
give an indication of what can be expected in more real- 
istic situations. 
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Karkowski et al., Fig. 1 



1.5e-17 



1e-17 



5e-18 



1e-18 



w=1.e7 f=0.2027 
w=5.e6 f=0.2045 
w=2.5e6 f=0.2031 




2.0001 2.0005 2.001 2.0015 2.002 2.0025 



(a) 



2.5 



2 



1.5 



0.5 

L 
4 

(b) 

600 r 

500 - 
400 - 
300 
200 
100 




w=0.0 f=1 .93e-4 
w=0.1 f=4.54e-4 
w=0.01 f=3.01e-4 



5 6 7 8 9 10 



w=1.e-2 f=3.71e-6 
w=1 .e-3 f=8.47e-6 
w=1.e-4 f=5.01 e-6 



20 40 60 80 100 

(c) 



20 



30 40 



200 300 400 



Fig.l. Initial energy densities of the electromagnetic field for various values 
of a ( (a): a=2; (b): a=4; (c): a=20) and w. The x-axis is in the LoglO - scale. 
The units of the y-axis are arbitrary. Numbers "f describe the ratio of the 
backscattered versus the initial energy. 
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Fig. 2. Initial energy densities of the gravitational field for various values of 
a ( (a): a=2; (b): a=4; (c): a=20) and w. The x-axis is in the LoglO - scale. 
The units of the y-axis are arbitrary. Numbers "f describe the ratio of the 
backscattered versus the initial energy. 



